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Properties of metals are fundamentally determined by their electron behavior, which is largely 
reflected by the electron work function (ip). Recent studies have demonstrated that many proper¬ 
ties of metallic materials are directly related to ip, which may provide a simple but fundamental 
parameter for material design. Since material properties are affected by temperature, in this ar¬ 
ticle a simple model is proposed to correlate the work function with temperature, expressed as 
<p(T) = ipo — 7 ^ , where 7 varies with the crystal structure. This tp-T relationship helps de¬ 
termine and understand the dependence of metal properties on temperature on a feasible electronic 
base. As a sample application, the established relationship is applied to determine the dependence 
of Young’s modulus of metals on temperature. The proposed relationship is consistent with experi¬ 
mental observations. 
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INTRODUCTION 


Material properties are fundamentally correlated to 
the electron behavior, which is largely reflected by the 
electron work function [1—8]. This correlation with work 
function includes a number of factors, including the 
Young’s modulus, thermal expansion, and heat capacity 
.1]. Recent studies [ —12] have shown that it may be 
more feasible to use the work function in material design 
compared to relevant quantum theories, since the latter 
are rather difficult to apply in material design, especially 
for structural materials. 

Many properties of materials are strongly affected by 
temperature. This is probably related to the influence of 
temperature on the behavior of electrons. The main ob¬ 
jective of this work is to establish a relationship between 
the work function and temperature, so that the depen¬ 
dence of the material properties on temperature can be 
predicted via the effect of temperature on the work func¬ 
tion, which also helps our fundamental understanding of 
such dependence. With the established ip — T relation¬ 
ship, we have predicted the dependence of Young’s mod¬ 
ulus on temperature as a sample application. 

Regarding the effect of temperature on Young’s modu¬ 
lus, a model to describe the variation in elastic modulus 
with temperature was first proposed by Born and Huang 
in 1954 [13], in which the Temperature dependence of 
elastic modulus results from non-harmonic changes in 
lattice potential energy, originally caused by lattice vi¬ 
brations [13, 14]. They demonstrated that the Young’s 
modulus varied with T 4 . Although consistent with the 
third law of thermodynamics, experimental results show 
that this dependence is quite limited as it is only valid 
when the temperature approaches absolute zero, which 


is ultimately negligible when considering larger ranges 
of temperature (e.g, 100 °K to 1000°A') that are more 
meaningful to engineering applications of metallic mate¬ 
rials at various temperatures. Wachtman et. al [15] have 
shown, through experimental data fitting, that Young’s 
modulus may be described as 


E = E 0 - BTexp(-T 0 /T) (1) 

where Eq, B and To are empirical constants. However, 
such an empirical equation does not provide a clear mech¬ 
anism for the described relationship. Besides, Anderson 
[ ] has shown that this equation is only valid when the 

variation of the Poisson ratio with temperature is small. 

A recent study by Hua and Li [9] has correlated the 
Young’s modulus of metals to their work function, with 
which the dependence of Young’s modulus on tempera¬ 
ture could be predicted if the effect of temperature on 
work function can be established. As demonstrated, the 
Young’s modulus of metals, E, has a sextic relation with 
the work function expressed as 


E 
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where a is Madelung’s constant, a is the equilibrium lat¬ 
tice constant, m is the electron mass, e is the elementary 
charge, and £q is the vacuum permittivity. A general 
expression that correlates the Young’s modulus with the 
work function is given as 


E = P (3) 

where ft = 0.02233^f is the average value for various 
crystal structures [ )]. This relationship is illustrated in 
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FIG. 1. Correlation between electron work function and 
Young’s modulus of metals (adapted from reference [9]). The 
collected experimental data points are around an average 
curve of E — [3<p 6 with /3 = 0.02233^rJ. The upper and 
lower bounds of the band, within which all the data points 
are located, correspond to /3 equal to 2/3 and 0.5/3, respec¬ 
tively. 


Fig.l for various metals. Alternative curves with differ¬ 
ent values of /3 for different crystal structures are demon¬ 
strated in the same figure. The region outlined by these 
curves represent the variation of the value for /3 from 0.5/3 
to 2/3. The value of j3 varies with the crystal structure. 

Thus, if the dependence of work function on tem¬ 
perature is established, the effect of temperature on 
Young’s modulus can be predicted. The work function- 
temperature relationship is of significance not only to 
Young’s modulus but also to other properties such as 
thermal expansion coefficient and heat capacity, etc. In 
addition to prediction of the dependence of intrinsic prop¬ 
erties of materials on temperature, ip(T) could also play 
an important role in acting as an alternative parameter 
for material design. 

The work function can be regarded as a barrier for the 
electrons to be moved from inside a solid to a point in 
vacuum immediately outside the solid surface [17]. In¬ 
creasing temperature should decrease this potential bar¬ 
rier, since electrons would be thermally excited and easier 
to be moved as the temperature rises. 

Dushman [18] has shown that the electron emission 
from metals is enhanced with an increase in tempera¬ 
ture. This is attributed to the energy absorption by 
electrons inside a metal, with an assumption that the 
increase in energy of each electron with temperature is 
equal to ^k B T (k B is the Boltzmann’s constant). Thus, 
the work function may be lowered as 


where n is the free electron density [ ]. The term | k B T 
is only the average energy of electrons according to the 
equipartition theorem, which states that the energy is 
shared between all accessible degrees of freedom, where 
each degree of freedom contributes }jk B T to the average 
internal energy [19]. Each degree of freedom is a param¬ 
eter that contributes to the state of a system and can 
be regarded as a classical harmonic oscillator with en¬ 
ergy k B T. However the equipartition theorem can only 
be used when the frequency of each oscillation is less 
than [IS]. According to Bardeen and Pines [ ], the 
plasma frequency for electrons (in SI units) in a solid is 
given by 


UJp — 



( 6 ) 


which is much larger than a t temperatures where 
materials are applied. Therefore the equipartition theo¬ 
rem is not applicable for the present case. 

According to the free electron model, the total energy 
of free electrons inside a metal is given by: 

/»00 

U(T) = / ef{e)D{e)de (7) 

Jo 

where f(e) is the Fermi-Dirac distribution, given by 


/(e) = 
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( 8 ) 


is the density of states and N is given as 


(9) 
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Equation (9) thus becomes: 


dN 3 N 
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( 11 ) 


The average energy for one electron at temperature T is 
therefore equal to 


<P = Vo - \ k sT (4) 

where, is the work function at T = 0 °K which is 
expressed as [9] 
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the integral in ( ) can be evaluated 
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The average values calculated using these two ap¬ 
proaches (i.e the equipartition theorem and the free elec¬ 
tron model) are not appropriate choices for determining 
the change in work function with temperature, since elec¬ 
trons with the highest energy, rather than the average 
energy, contribute to the change in work function. 


DERIVATION OF ip(T) USING LENNARD JONES 
POTENTIAL 

The potential between two atoms, according to the 
Lennard-Jones potential is commonly expressed as [10] 

V(r)=s 6 [(^) 12 -2(^) 6 ] (14) 

where £5 is the maximum value for the depth of the po¬ 
tential well, r is the distance between the two atoms and 
r e is the equilibrium distance. It has been shown that 
the electron work function is related to the bond energy 
[10]. This relationship is given as 

tp(r e ) = Cs 1 / 6 (15) 

where v( r e) is the work function at equilibrium and C is 

a constant of proportionality. Since V(r e ) = — 

V{r e ) = -^£ 5 b /6 ip(r e ) (16) 

Because the atomic displacement or vibration, Ar = r — 

r e , is within a small range (r < 1.1 r e , otherwise the 
bond would become unstable [ 21 ]), equation (16) can be 
written as 


V(r) = -^£)i /f V(r) (17) 

Combining equations (14) and (17), the expression for 
the work function is finalized as: 

v( r ) = ~~v( r e)[( — ) 12 ~ 2 ( —) 6 ] (18) 

r r 

Letting x = <r ~ re> , v( r ) can be written as 


x = where we may have g = 252^ and 

/ = 36Thus, x becomes 


x= ( 


7 k B T 
48 ip 0 
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Combining equations (19) and (20), we have 


( 20 ) 


<p(T)=<Po~ 36(^^f— )Vo ( 21 ) 

4 » ip 0 

This expression can be generalized for solids by consider¬ 
ing the potential due to the interaction of other adjoining 
atoms. The second term in equation ( 9) after the neg¬ 
ative sign can be generally described as: 
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= 7^o( — ) 2 (22) 

Vo 


where < Ti — r ei > is average deviation or displace¬ 
ment of atom i away from the equilibrium position (cho¬ 
sen as origin). w t is an energy contribution factor that 
counts the influence of the zth neighbor atom on the 
square of relative oscillation amplitude, x 2 . This con¬ 
tribution is distance-dependent, e.g., the contributions 
from the nearest neighbour, 2 nd and 3 rd neighbour atoms 
are w r = 1 , w r = 3.1 x 10 “ 2 and w r = 2.7 x 10 —3 , 
respectively (determined based on the Lennard-Jones po¬ 
tentials for r e , 2 r e , 3 r e , ... ). <p(T) then becomes 


(kBT ) 2 

V(T) = Vo~ 7- (23) 

Vo 

Combining equations (3) and (23), the Young’s modulus 
becomes 


E = /%o - 7 ^^] 6 (24) 

Vo 

Since the values for work function are available at the 
room temperature (i.e 295 K), the relationship for the 
Young’s modulus can be adjusted as 


v(r) = ~Vo[(l + x) 12 — 2(l + a;) 6 ] « <p 0 -36a;Vo (19) 

Since the distance between the two atoms (the bond 
length) is affected by temperature, an expression for a 
temperature-dependent work function can be derived us¬ 
ing this relationship. Since the change in potential en¬ 
ergy is similar to that in work function (see Eqs. (14) 
and 0 8)), the temperature dependent x varies as [ 


E = p[v 295 - 7 


(kBT ) 2 

V295 


(25) 


where r = T — 295A7 

7 in eq. ( 22 ) can be written as 7 = £( jg ) 2 and £ is: 


z = 


36 

X? 


N 

y^ < r i ~~ r ei > )2 


(26) 
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Metal 

7 

p\p\ 

A1 

583 

0.5 

Fe 

438 

1.1 

Ag 

478 

0.6 

Ni 

318 

0.5 

Cu 

307 

0.6 


TABLE I. Calculated values for 7 and estimated values of /3, 
from figure 1 . 

where Yi = —-— TV| > = a^T and is the linear ther- 

r ei 

mal expansion coefficient of the metal. Considering the 
periodicity of the lattice, we have: 

r ei = r ei + (i - l)a (27) 

where a is the lattice parameter and r ei is the nearest 
neighbour distance, which equals ^ for face-centred cu¬ 
bic (fee) crystals and for body-centred cubic (bcc) 
crystals, respectively. Thus, £ for fee structures becomes: 


Iron Silver 



Temperature [°K] 

Nickel Copper 



FIG. 2. Variation in Young’s modulus of metals with temper¬ 
ature. The solid line represents the prediction and the dashed 
lines are the experimental values. 


£fcc 


36 \ - 


XI + - !) 

1 + y/2(i - 1) 
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Similarly, for bcc structures, we have 


(28) 


36 f Xi + ^(i-l), 

(ecc = ^( l + ^ {i _ iy 


(29) 


Considering the interaction between two neighbouring 
atoms and calculating xi a t the melting point (the dis¬ 
tance between two neighbouring atoms in a metal is max¬ 
imum at the melting point), the values for 7 for different 
metals are calculated and presented in Table 1. With 
the 7 values and those for /?, the variation of the Young’s 
modulus with temperature is plotted and presented in 
Fig. 2. Experimental data are cited from ref. [22-36]. 


SUMMARY AND DISCUSSION 

In summary, we propose a relationship between the 
work function of metals and temperature, ip(T ) = ipo — 
7 B ’ , where the coefficient 7 is dependent on the crys¬ 
tal structure. Based on this relationship, the temperature 
dependence of Young’s modulus is established. Using 
iron, silver, aluminum, nickel and copper as examples, 
variations in their Young’s moduli with temperature were 
predicted, which are supported by reported experimen¬ 
tal observations. The <p-T relationship is of significance 
not only to Young’s modulus but also helps predict the 


dependence of other intrinsic properties of metals on tem¬ 
perature on a feasible electronic base. The proposed rela¬ 
tionship is general, since there is no specific assumption, 
which could limit the applicability of the relationship, 
required in the derivation. 
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